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In the present study, we performed molecular-dynamics simulations and investigated dynamical heterogeneity 
in a supercooled liquid under a steady shear flow. Dynamical heterogeneity can be characterized by three 
quantities: the correlation length S,i{t), the intensity Xi{t), and the lifetime Thetero(0- We quantified all 
three quantities by means of the correlation functions of the particle dynamics, i.e., the four-point correlation 
functions, which are extended to the sheared condition. Here, to define the local dynamics, we used two 
time intervals t = Ta and T,igp; Tq, is the a- relaxation time, and T,igp is the time at which the non-Gaussian 
parameter of the Van Hove self-correlation function is maximized. We discovered that all three quantities 
(?4(0: Xi{i)^ and Thctcro(O) dccrcasc as the shear rate 7 of the steady shear fiow increases. For the time 



interval t — Ta, the scalings ^4(tq,) 



7 



-0.08 



-0.26 



and Thctcro{Ta) ~ 7 were obtained. The 



steady shear fiow suppresses the heterogeneous structure as well as the lifetime of the dynamical heterogeneity. 
In addition, our results demonstrated that the a-relaxation time Tq dependences of three quantities coincide 
with those at equilibrium. This means that all three quantities of the dynamical heterogeneity can be mapped 
onto those in the equilibrium state through Tq,. 



PACS numbers: 64.70.P-, 61.20.Lc, 83.50.Ax 



I. INTRODUCTION 

When liquids are cooled toward the glass transi- 
tion temperature Tg, the dynamics not only drasti- 
cally slow down, but also become progressively more 
heterogeneous^"^. "Dynamical heterogeneity" has at- 
tracted much attention because it may lead to a bet- 
ter understanding of slow dynamics^"'^^. Dynamical het- 
erogeneity is characterized by three quantities: the cor- 
relation length ^4, the intensity X4, and the lifetime 
Thotoro^^"^^- The correlation length ^4 and the intensity 
X4 characterize the static properties of the dynamical het- 
erogeneity and can be quantified by means of the static 
structure factor of the particle dynamics, i.e., the so- 
called four-point correlation function^^"^^. The lifetime 
Thetero measures the dynamical properties of the dynam- 
ical heterogeneity, and we can evaluate the lifetime by 
the time decay of the correlation function of the parti- 
cle dynamics^^"^^. In our previous studies^^'^"^, we have 
consistently quantified these three quantities (correlation 
length, intensity, and lifetime) from a single order param- 
eter and its correlation functions. The length scale and 
the time scale of the dynamical heterogeneity tend to di- 
verge as the temperature approaches the glass transition 
point. 

Dynamical heterogeneity has also been investigated for 
a sheared situation^*"^^. When applying a steady shear 
flow to supercooled liquids, marked shear-thinning can 
be observed^'^"'^^ . In the sheared situation, the dynamics 
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become not only much faster but also more homogeneous; 
therefore, the shear fiow suppresses the dynamical hetero- 
geneity. Refs.^^"^^ found that the intensity X4 of the dy- 
namical heterogeneity decreases with an increasing shear 
rate 7 as X4 ~ 7"'" with ^ = O.S^^'^" or ^ = 0.4 - 0.6^^ 
To the best of our knowledge, the lifetime Thctcro of the 
dynamical heterogeneity has not yet been examined for 
the sheared situation. Dynamical heterogeneity may play 
an important role in the drastic change of the dynamics 
due to not only the decreasing temperature but also the 
increasing shear rate. However, the role of dynamical 
heterogeneity remains unclear and is an important ques- 
tion. 

Even in a strongly sheared state, supercooled liq- 
uids exhibit an almost isotropic structure and dy- 
namics that can be captured via two-point correlation 
functions^*'^^''^^. This is in marked contrast to common 
complex fiuids, such as polymer solutions and colloidal 
suspensions, which exhibit anisotropic dynamics or struc- 
tural changes that are induced by the shear fiow'^^"'^*. 
The isotropic feature of supercooled liquids strongly sup- 
ports the possibility of the "mapping concept" , that is, 
that the sheared non-equilibrium state can be mapped 
onto the quiescent equilibrium state. In fact, previous 
studies^^'^^ have demonstrated that the shear viscosity 
77 and the diffusion constant D in the sheared state can 
be mapped onto those in the equilibrium state through 
the a-relaxation time Tq. The research by Ref.'"' also 
reported that the shear stress and the inherent structure 
energy can be mapped through the effective tempera- 
ture Toffj which is measured from the relation between 
the static linear response and the variance of the pres- 
sure. More interestingly, there are some results regard- 
ing the mapping of the dynamical heterogeneity. An ear- 
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lier study^ examined the heterogeneous structure of the 
bond breakages under the sheared situation and found 
that tlic correlation length of the bond breakages is de- 
scribed as a function of only the a-relaxation time. Fur- 
thermore, recent research^" has found that the intensity 
X4 can be scaled as n* ^ (-yAc/)^)""-^, where 7 is the 
shear rate of the shear flow, and is the distance from 
the random close-packing volume fraction. This result 
means that the intensity can be mapped through a sin- 
gle parameter •yAcj)'^. 

In the present study, we examined the dynamical het- 
erogeneity in the sheared non-equilibrium situation using 
molecular-dynamics (MD) simulations. Three quantities 
(correlation length, intensity, and lifetime) characterizing 
the dynamical heterogeneity were quantified by means of 
the four-point correlation functions. We found the be- 
haviors of decreasing of not only the correlation length 
^4 and the intensity Xi but also the lifetime Thctcro with 
increasing the shear rate 7. Furthermore, we also ex- 
amined the validity of the mapping concept for all three 
quantities of the dynamical heterogeneity. 

The mapping concept can lead to the concept of an ef- 
fective temperature, which is also interesting and impor- 
tant for glassy systems. For glassy systems, it has been 
suggested that when considering the long time scale, the 
equilibrium form of the; FDT holds well with the temper- 
ature T, which is rcplaccid by a different value denoted 
as Teff, i.e., the effective temperature**^"^^. If the concept 
of an effective temperature is valid, then the effective 
temperature plays a role as the temperature in the non- 
equilibrium situation, and the sheared non- equilibrium 
state can correspond to the equilibrium state with the 
effective temperature Tcff- However, there are difficulties 
regarding the mapping concept and the effective temper- 
ature. A recent numerical simulation^"^ found out the 
anisotropic dynamics in a sheared supercooled liquid by 
the four-point correlation function. This finding indicates 
that the sheared state cannot completely be mapped onto 
the equilibrium state. Furthermore, the validity of the 
value Teff, which is calculated from the relation between 
the response fimction and its associated correlation func- 
tion, remains uncertain. Although several different ob- 
servables yield a common value of Teff^^~*^'^®, some ob- 
servables do not provide consistent values of Tcg^^'^^. 
Therefore, only one value of Tea may not completely de- 
scribe the non-equilibrium state. Despite these difficul- 
ties, the mapping concept and the effective temperature 
are very attractive and interesting. 

The paper is organized as follows. In Sec. II, we briefly 
review our MD simulation and the correlation fun(;tions 
of the particle dynamics. In this section, we extend the 
correlation functions to the present sheared condition. In 
Sec. Ill, we present our results. First, we show the shear 
rate dependences of the correlation length, the intensity, 
and the lifetime, and then we demonstrate whether the 
mapping concept is valid for these three quantities . In 
Sec. IV, we summarize our results. 



II. SIMULATION MODEL AND CORRELATION 
FUNCTIONS OF PARTICLE DYNAMICS 

A. Simulation model 

We performed MD simulations in three dimensions. 

Our simulation model is a mixture of two different size 
atomic species, 1 and 2. The particles interact via the 
soft-sphere potential Vab{r) = e{aab/rY'^ with Uab = 
{da + (Jb)l'2,^ where r is the distance between two par- 
ticles, Ua is the particle size, and a,b Cz 1,2. The interac- 
tion was truncated at r = 3aab- We took the mass ratio 
to be 1712 /mi = 2 and the size ratio to be a2/<Ji = 1.2. 
This diameter ratio avoided system crystallization and 
ensured that an amorphous supercooled state occurred 
at low temperatures^^. As in our previous study we 
used two systems: a small system with Ni = N2 = 5x 10^ 
(N = Ni_ + N2 = 10'*) particles and a large system with 
TVi = iV2 = 5 X 10'' {N = Ni + N2 ^ 10^') particles. 
We quantified the correlation length ^4 (t) and the inten- 
sity X4(t) with the large system, and we quantified the 
lifetime Thctcro (i) with the small system. In the follow- 
ing, space, time, and temperature were measured by tJi, 
To = (micrf /e)^/^, and e/fcs, respectively. The particle 
density was fixed at a value of p = (TVi -I- N2)/V = 0.8, 
and the system length was L = V^l"^ = 23.2 and 50.0 for 
the small and large systems, respectively. The temper- 
ature was set as T = 0.772, 0.473, 0.352, 0.306, 0.267, 
and 0.253. Note that the freezing point of the correspond- 
ing one-component model is approximately T = 0.772^"'^. 
At r = 0.253, the system is in a highly supercooled state. 
We applied a steady shear flow on our system with the 
Lees- Edwards periodic boundary condition^^. We inte- 
grated the so-called SLLOD equations of motion with 
the Lees-Edwards periodic boundary condition, and the 
temperature was maintained at a desired value using a 
Gaussian constraint thermostat. Here, we set the x axis 
and the y axis along the flow direction and the velocity 
gradient direction of the steady shear flow, respectively. 
The mean velocity profile is (v) = jye^ , where is the 
unit vector in the x direction. The details of this simu- 
lation model can be found in previous studies^^'^^. 

B. Correlation functions of particle dynamics 

In our previous studies^®'^^, we investigated dynam- 
ical hetc;rogc>iieity in the unsheared equilibrium state. 
With the correlation functions of the particle dynamics, 
we quantified the correlation length £,4{t), the intensity 
X4(i), and the lifetime Thctcro (i), three of which char- 
acterize the static and dynamic properties of dynami- 
cal heterogeneity. In those studies^^'^'', we defined an 
order parameter SDi{r,to,t) and its Fourier component 
d'Di{q,to,t) as the local fluctuations in the particle mo- 
bility. In the present study, we extended 6'Di{r, to, t) and 
5'Di{q,to,t) to the sheared condition as follows^". A po- 
sition r{to) = {x{to),y{to)) on the reference frame at a 
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FIG. 1. Schematic illustration of the time configuration of the 
correlation functions of the particle dynamics: (a) Si{q,t) (b) 
S'D{q,ts,t). 



time to moves to T(to) +ity{tQ)ex at a time tQ + thy the 
mean shear flow. Therefore, we must define the mobility 



alAto,t) of the particle j by 



0-ij[to,t) - — — (ij 



{[rijita + t) - rij(to) - jtyij{to)e^]^) ' 



Note that when alj{to, t) > 1 {alj{to, t) < 1), the particle 
j moves more (less) than the mean value of the single- 
particle displacement, i.e., the particle j is mobile (im- 
mobile). Using this definition of particle mobility, we ex- 
tended 5T>i{r,to,t) and 5'Di{q,to,t) to the sheared con- 
dition as 

SA{r,to,t)=Y,[alj{to,t) - l]Sir~Ri^{to,t)), 
j=i 

Ni 

6T>i{q,to,t) = ^[alj{to,t) - 1] exp[-zq • Rij{tQ,t)], 

3 = i 

(2) 

where Rijito,t) = {rij{to) + rij{to + t))/2. 

Using the order parameter d'Di{q,to,t), we define the 
spatial correlation function S4{q,t) of the particle dy- 
namics as 

S^iq^t) = ■^{SVi{q,0,t)SV,i-q,0,t)), (3) 
iVi 

where the (• • • ) represents the ensemble average over time 
and the angular components of the wave vector q. 
The function S4^{q,t) represents the spatial correlation 
of the particle dynamics in the time interval [0,t\. The 
time configuration of S'4(g, i) is schematically illustrated 
in Fig. 1(a). We are able to quantify the correlation 
length £,4{t) and the intensity Xii^) by fitting S4{q,t) to 
the Ornstein-Zernike (OZ) form at small wavenumbers 



Si{q,t) 



(4) 



Note that to obtain accurate values of £,i{t) and X4(*)j 
Si{q,t) must be fitted to the OZ form in the range of 
qS,i{t) < 1.5^°'^^. We used a large system with 10^ par- 
ticles and carefully fitted Si{q,t) to the OZ form. 

In addition, we also define the time correlation function 
Svil^tsji) of the particle dynamics as 

Sviq,ts,t) = {SVi{q,ts +t,t)SVi{-q,0,t)), (5) 

where the (• • • ) represents the ensemble average over the 
initial time and the angular components of the wave 
vector q. The function Sviqjtsjt) represents the corre- 
lation of the particle dynamics between two time inter- 
vals [0,t] and [ts + t,ts + 2t\, where the value tg is the 
time separation. The time configuration of 6*15(5, is, t) is 
schematically illustrated in Fig. 1(b). Here, when calcu- 
lating the ensemble average over the wave vector q, we 
took only q perpendicular to x direction, i.e., q perpen- 
dicular to the shear fiow direction, to eliminate the time 
decay due to the shear flow. As the time separation tg 
increases, S'D{q,ts,t) decays in the stretched exponential 
form, 



ST>iq,ts,t) 
Sv{q,0,t) 



exp 



Th{q,t) 



(6) 



where Th{q,t) is the wavenumber-dependent relaxation 
time of Sx>(q,ts,t). The lifetime Thetero(i) was deter- 
mined as Thotoro(^) = Th{q^t) aX q = 0.38. We used a 
small system with 10"* particles to obtain the lifetime 

'^hctcro(^) ■ 

As we mentioned previously, recent research^'' has in- 
dicated that the correlation functions of the particle dy- 
namics (the four-point correlation functions) depend on 
the direction of the wave vector q in the sheared situa- 
tion, although the two point correlation functions hardly 
exhibit any anisotropy. Therefore, three quantities (^4(^)1 
X4(t), and Thotoro(O) 9'lso depend on the direction of q. 
In the present study, we did not examine the directional 
dependence of the dynamical heterogeneity, even though 
it is interesting and important. Rather, we examined 
the dynamical heterogeneity and its three quantities av- 
eraged over the direction of the wave vector q. 



III. RESULTS 

To define the particle dynamics, we used two time in- 
tervals as in our previous studies^^'^^. One is the a- 
relaxation time Tq, defined by Fs{km,Ta) — e^^. Here, 
Fs{k,t) is the self-part of the density correlation func- 
tion for Particle Species 1, which is defined as Fs{k,t) — 

lM(Efiie"''''''''^^*^"'"'^'^°^~^*'''^^*^*'^')> and k„^ = 
27r is the first peak wave number of the static structure 
factor. We note that, as already shown in Refs.^*''^^'^^, 
any angular dependence in Fs{k,t) is barely notable de- 
spite the presence of the shear flow in the x direction. 
Figure 2(a) shows Tq, as a function of the shear rate 7. 



10"^ 10'^ 10"* 10'^ 10"^ 10''' lo" 

7 

FIG. 2. The shear rate dependence of (a) Ta and (b) Tngp for 
various temperatures. We used these two time intervals to 
define the particle dynamics. 



The value monotonically decreases with increasing 7 
as Ta ~ 'j"" with f ~ 0.8^*'^^'^^. Note that the marked 
shear thinning occurs in the region of the shear rate that 
we considered. 

The other time interval is the time r^gp at which 
non-Gaussian parameter 012 (i) of the Van Hove self- 
correlation function is maximized. Here, we introduce 
a displacement Ar[j{t) of the particle j, which is defined 

as Ar[^{t) = rij{t) - rij(O) - 7 J* dsyij{s)e^, in which 
the contribution from the convective transport by the av- 
erage shear flow is subtracted^^. By using this displace- 
ment, we define 03 (t) = 3{[Ar[{t)]^) /5{[Ar[{t)]Y - 1- 
In Fig. 2(b), we plot Tngp as a function of 7. Similar to 
'Ta, Tngp also dccrcascs with increasing 7 as Tngp ~ ^-o.65^ 
This behavior of r„gp was also experimentally observed 
in Ref.^^ 



A. Heterogeneous dynamics under shear flow 

First, we visualized the heterogeneous dynamics under 
a steady shear flow. We show the spatial distribution of 
the particle mobility a1j{to,t) for Particle Species 1 in 
Fig. 3, where the particles are drawn as spheres with 
radii a^{to,t) located at Rj{to,t). The time interval is 
t ■= Ta and t = Tngp, and the temperature is T = 0.267. 



(c) t = Tngp and 7 = 10 ^ (d) < = Tngp and 7 = 10 



FIG. 3. The visualization of the heterogeneous dynamics for 
Particle Species 1. The temperature is 0.267. The time inter- 
val is [to,to-\-Ta\{t = Ta) iu (a) aud (b) and [to, io + Tngp](t = 
Tngp) in (c) and (d). The shear rates are 10~^ in (a) and (c) 
and 10~^ in (b) and (d). The radii of the spheres are a\j{to,t), 
and the centers are at Rij {to , t) . The red and blue spheres 
represent a\j{to,t) > 1 (mobile particles) and aij{to,t) < 1 
(immobile particles), respectively. 

We demonstrate that the heterogeneity is significant at 
7 = 10~^ but greatly weakened at 7 = 10""'^ for both 
time intervals t = Ta and Tngp. Therefore, the strong 
shear flow suppresses the heterogeneous structure. The 
same observation has also been observed in the hetero- 
geneity of the bond breakage^® and of the "overlapping" 
function^^"^^''^°. 



B. Correlation length and intensity of dynamical 
heterogeneity 

Next, we quantified the correlation length ^4{t) and the 
intensity X4{t) of the dynamical heterogeneity. We cal- 
culated the spatial correlation function 5*4(5, t) defined 
in Eq. (3). Note that we used a large system with 10^ 
particles for the calculation of S4{q,t). Figure 4 shows 
the wave number q dependence of Si{q,t) for various 
shear rates 7. The temperature is T = 0.267. It can 
be observed that S4{q,t) at small wave numbers q (long- 
distance scales) for t = Ta and Tngp decreases as the shear 
rate 7 increases. Therefore, the heterogeneous structure 
is weakened due to the steady shear flow as is also shown 
in Fig. 3. 

We quantified the correlation length £,4{t) and the in- 
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FIG. 4. The spatial correlation function S4{q,t) for Parti- 
cle Species 1. The temperature is T = 0.267. The time 
interval t is Ta in (a) and Tngp in (b). The shear rate 7 is 

0, 10~'^ 10^", 10"'', 10"^ and 10"^ from the highest curve 
to the lowest. Note that 5*4(9, i) was calculated with a larger 
system N = . 



tensity X4(t) by fitting the function S4{q,t) to the OZ 
form in Eq. (4). Note that S4{q,t) was carefully fitted 
to the OZ form in the range of (7^4 (t) < 1.5 to obtain 
accurate values of S,4{t) and X4(i)^"'^^- In Figs. 5 and 6, 
we show the shear rate 7 dependences of the correlation 
length ^4{t) and the intensity Xiit) for various temper- 
atures. It is clear that both ^4{t) and X4{t) decrease 
as 7 increases for the time intervals t = and Tngp. 
For the time interval t = t^, ^4(rQ,) and Xi{Ta) can be 
scaled as S,4(Ta) ^ 7^0-"8 ^j^j X4(tq,) ^ ^^-0.26 pjgs_ 
5(a) and 6(a), respectively. The scaling component 2.6 
of X4('ra) ^ 7"°-26 agrees well with 0.3 in Refs.^^'^° but 
disagrees with 0.4 — 0.6 in Ref.^^. 



C. Lifetime of dynamical heterogeneity 

Next, we quantified the lifetime Thetero(i) of the dy- 
namical heterogeneity in terms of the correlation function 
S-D{q,ts,t) defined in Eq. (5). We calculated ST>{q,ts,t) 
for the time intervals t = Ta and Tngp. Figure 7 shows 
the time decay of S-piq, ts,t) a.t q = 0.38 for various shear 
rates 7. Note that q = 0.38 is the smallest waver num- 



FIG. 5. The shear rate dependence of the correlation length 
^4{t) for various temperatures. The time interval t is (a) Ta 
and (b) Tngp. 
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FIG. 6. The shear rate dependence of the intensity x^it) for 
various temperatures. The time interval t is (a) Ta and (b) 

Tngp. 
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FIG. 7. The time decay of ST){q,ts,t) for Particle Species 
1 at (/ = 0.38. The time interval t is Tc in (a) and Tngp 
in (b). The temperature is 0.267. The shear rates 7 are 
0, 10"^, IQ-*, 10"^, 10-^ and 10"^ from right to left. 
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FIG. 8. The shear rate dependence of the lifetime rhctoro(t) 
for various temperatures. The time interval t is (a) Ta and 

(b) Tngp. 



bcr (the longest distance scale) that we used in calculat- 
ing Sxi{q, ts,t). Sxi{q_, ts,t) decays in the stretched expo- 
nential form with increasing the time separation tg- As 
shown in Fig. 7, as the shear rate 7 increases, /St? (9, ts,t) 
for both t = Ta and Tngp decays faster. 

Wc quantified the lifetime Tiictcro(t) using the relax- 
ation time of Sx>{q,tg,t) as in Eq. (6). Figure 8 shows 

'^hcicro 

{t) for t = Ta and Tngp as a function of 7. Similar 
to Ta and Tngp, Thctcro(i) also becomes small as 7 be- 
comes large. Here, the scalings Thetero(T"a) ~ -y~°-^^ and 
Thetero(Tngp) ~ 7~°'^* are observcd. The steady shear 
flow decreases the lifetime of the dynamical heterogene- 
ity as well as the correlation length and the intensity. By 
comparing three scalings ^4 ~ -y"'^-'^^, X4 ~ t'"^-^^, and 
Thetero ^ ^"O.ss ^^iQ time interval t = Ta, we note 
that the lifetime is much more sensitive to the shear rate 
7 than the correlation length and the intensity. Thus, 
the dynamical properties of the dynamical heterogeneity 
are altered much more by the shear flow than the static 
properties of the dynamical heterogeneity. 



D. Mapping concept for correlation length, intensity, and 
lifetime of dynamical heterogeneity 

Finally, we demonstrated that three quantities, the 
correlation length ^i{t), the intensity X4(^)) and the life- 
time Thetero(i), cau be mapped onto those in the equi- 
librium state. In Fig. 9(a), we plot £,4{t) versus Tq, for 
both time intervals t = Ta and Tngp. In the same fig- 
ure, we also plot the values in the equilibrium state. It 
can be observed that the values ^i{t) for t = Ta and 



'ngp 



collapse onto a single curve as a function of Ta- 
This result means that the correlation length ^i{t) in 

the sheared state can be mapped onto that in the equi- 
librium state through the a- relaxation time Tq. For the 
time interval t = Ta, we obtained the scaling relation- 
ships between ^4{Ta) and Tq,: ^4{Ta) ~ t°'-^ (a power law 
fit) or ^4(tq) ^ (lnTQ,)°'^° (a slower-than-power law fit). 
These fits are exactly the same as the equilibrium fits 
within the numerical errors'^^. In addition, in Fig. 9(b), 
we show X4:{t) versus ^4{t) for t = Ta and Tngp. From 
Fig. 9(b), similarly to £,4{t), Xiit) a-lso collapses onto a 
single curve as function of ^4{t). Therefore, the values 
Xiit) in the sheared state can also be mapped onto those 
in the equilibrium state. For the time interval t = Ta, 
we obtained the scaling relationship between X4(Ta) and 
^4(Ta). It is a power law fit X4(tq) ~ ^4(^0)^'^, which is 
same as the equilibrium fit within the numerical errors^''. 
It should be noted that Ta decreases with increasing 7 
as Ta - 7""■^ which leads to ^4(ra) - t^-^ - 7-"-"*^ 



and X4(ra) ~ U{Ta) 



3.3 



1 



-0.26 



These scaling relations 



coincide with the observations in Figs. 5(a) and 6(a), 
respectively. 

Furthermore, figure 10 shows the lifetime Thetero(i) 
versus Ta for t = Ta and Tngp. This figure shows 
that Thetero(Ta) and Thctcro (^ngp ) fall into the lines 

Thetero(Ta) ~ T^'^° and Thetero(Tngp) ~ T^^^, similarly tO 
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FIG. 9. (a) The correlation length ^4(f) versus the a- 
rclaxation time Tc. (b) The intensity X/^(t) versus the cor- 
relation length ^4(t). The time intervals are t = Ta and Tngp. 
The white symbols arc> tlie values of ^4(i) and X^i't) under 
shear flows. The black symbols are the values in the equilib- 
rium states from Ref.^^. The straight lines are power law fits: 
e4(r«) ~ r«-i±0-«i in (a) and X4(rc) ~ ^rcf'^^^-^ in (b)- 
The dashed curve is fit to Uiro,) ~ (InTc,)"-**'' "^ 
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FIG. 10. The Ufetime (t) versus the a-relaxation time 

Ta- The time intervals are t = Ta and Tngp. The white sym- 
bols are the values of rhotcro(t) under shear flows. The black 

symbols are the values in the equilibrium states from Ref.^^ 
The straight lines are power law fits: Thetero(Ta) 

„v,^ „ ^ ,0.92±0.02 

and Theterol.T-ngpj ~ Ta 



^1.10±0.02 



^4 (t) and X4 (t) . These lines are exactly same as the equi- 
librium ones within the numerical errors^^. We note that 



^1.10 



-0.88 



and Thctcro(T, 



ngp; 



^0.92 



Thctcro 

^-0.74 coincide with the observation in Fig. 8. There- 
fore, not only the correlation length and the intensity 
but also the lifetime can be mapped onto those in the 
equilibrium state through the a-relaxation time. As we 
mentioned, previous studies^^'^^'*° have found that the 
material properties, such as the shear viscosity rj and 
the diffusion constant D, under steady shear flow can be 
a single function of the a-relaxation time or the ef- 
fective temperature Toff. Our results demonstrated that 
the mapping concept is valid for the material properties 
(the shear viscosity and the diffusion constant) as well 
as the properties of dynamical heterogeneity (the corre- 
lation length, the intensity, and the lifetime). 



IV. CONCLUSION 

In conclusion, we examined the dynamical heterogene- 
ity in a supercooled liquid under sheared conditions. 
Dynamical heterogeneity can be characterized by three 
quantities: the correlation length ^4(t), the intensity 
X4:{t), and the lifetime Thotoro(0 • Using the correla- 
tion functions of the particle dynamics, we quantified all 
three quantities for two time intervals t = Ta and Tngp. 
Our results demonstrated that not only the correlation 
length and the intensity but also the lifetime decrease as 
the shear rate j increases. Therefore, the steady shear 
flow suppresses the length scale as well as the time scale 
of the dynamical heterogeneity. For the time interval 
t = Ta, three quantities can be scaled as ^4 ~ j-o.08^ 
X4 ^ -y"*^-^^, and Thetero ~ j-o-ss^ which indicates that, 
due to the steady shear flow, the lifetime decreases much 
more than the correlation length and the intensity. 

Furthermore, we numerically demonstrated that the 
values of ^4{t), X4(f), and Theteio{t) collapse onto a sin- 
gle curve as function of the a-relaxation similarly to 
the material properties such as the viscosity rj and the 
diffusion constant D. For the time interval t = Ta, we 
obtained exactly the same scaling relations between the 
three quantities as the equilibrium ones, i.e., ^4 ~ t^'^ 
(or ^ - (lnr„)0-8°), X4 ~ ^|■^ and Thetero ~ ri io. Our 
results indicate that three quantities of the dynamical 
heterogeneity can be mapped onto those in the equilib- 
rium state, i.e., the mapping concept for three quantities 
is valid. 

In the present study, we did not acknowledge the 
anisotropy of the dynamical heterogeneity, and the values 
of ^4(^)1 X4{t), and Thetero (i) wcre averaged over the di- 
rection of the wave vector q. However, as we mentioned, 
the dynamical heterogeneity becomes anisotropic in the 
sheared situation, and the anisotropic dynamics can be 
related to shear-thinning^°. How the dynamical hetero- 
geneity and its three quantities depend on the direction 
of the wave vector q is interesting and important subject. 
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